INFINITE ENERGY SOLUTIONS FOR DAMPED NAVIER-STOKES 

EQUATIONS IN M^ 
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SERGEY ZELIK 

Abstract. We study the so-called damped Navier-Stokes equations in the whole 2D space. 
The global well-posedness, dissipativity and further regularity of weak solutions of this problem 
in the uniformly-local spaces are verified based on the further development of the weighted 
energy theory for the Navier-Stokes type problems. Note that any divergent free vector field 
uo € L°°(]R'^) is allowed and no assumptions on the spatial decay of solutions as |2;| -^ oo are 
posed. In addition, applying the developed theory to the case of the classical Navier-Stokes 
problem in R^, we show that the properly defined weak solution can grow at most polynomially 
(as a quintic polynomial) as time goes to infinity. 
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1. Introduction 
We study the following damped Navier-Stokes system in the whole space a; G M^: 



fl.l) 



dtu + {u, Vx)u 
div n = 0, 



t=o 






au + VxP + 9, 



where a is a positive parameter. These equations describe, for instance, a 2-dimensional viscous 
liquid moving on a rough surface and are used in geophysical models for large-scale processes 
in atmosphere and ocean. The term au paramejtemes the extra dissipation occurring in the 
planetary boundary layer (see, e.g., |T^J; see also fSr^or the alternative source of damped Euler 
equations). 

The mathematical theory of d|imjicd^N^ and Euler equations is of a big current 

interest, see l^"iS7iS7iinriII*"LEj and references therein. However, most part of these papers study 
either the case of bounded underlying domain (e.g., with periodic boundary conditions) or the 
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case of finite energy solutions in the whole space M? and very few is known about the infinite- 
energy solutions (for instance, starting with uq E L°°(M^) or from the proper uniformly local 
Sobolev space) which are natural from the physical point of view (for instance, for anphcations 
to the so-called uniform turbulence theory or/and statistical hydrodynamics, see \'lt\ for the 
details). Indeed, the main technical tool for the mathematical study of such solutions is the 
vorticity equation 

(1.2) dtto + {u,Vx)uj = AxU) — auj + cnvlg, a; = curln 

which possesses (in 2D) the maximum principle and, as a consequence, the dissipative L°°- 
estimate for the vorticity is almost immediate: 

(1.3) II curln(t)||Loo < || curluo||L°° H — 1| curlf^Hioo. 

a 

Nevertheless, estimating the L°° or uniformly bounded norms of the velocity field u based on the 
vorticity estimate (pi.3J) is far from being straightforward (in particular, due to the ill-posedness 
of the Helmholtz projectors in the L°°-type spaces) and, to the best of our knowledge, only 
exponentially growing in time estimates of u{t) can be obtained based solely on the vorticity 
estimate (pl.aj) . The situation becomes even worse in the case of classical Navier-Stokes problem 
(a = 0), where the vorticity estimate may grow linearly in time, but the best known estimate 
of the L°°-norm of u{t) is super-exponential: 

(1.4) ||n(t)|Uoo<Cie^2t2, 

|GMS01.ST07 IAM05 

where Ci and C2 depend on uq and g, see |n, iYj (see also |Tf7or the polynomial estimates in 
the case of periodic boundary conditions in one space direction). 

Thus, the alternative/complementary to the vorticity estimate ideas/methods are necessary 
to obtain more realistic upper bounds for the solution u{t). The most natural would be to use 
the so-called weighted energy estimates which are known as very effective tools in the general 
theory of dissipative PDEs in unbounded domains, see |T5,~'25'J^and the references therein. The 
main obstacle to the use of such estimates for the Navier-Stokes type problems is the fact that 
the inertial term {u,'Vx)u does not vanish being multiplied by the weighted term u<p (after the 
integration by parts, here ^{x) is a weight function), but produces the non-sign definite cubic 
term of the form ip'u"^. Since the other terms in the weighted estimate are at most quadratic, it 
was not clear how to control this higher order non sign-definite term and to close the estimate. 
Another problem is that wp is not divergence free, so the pressure term also does not vanish 
and requires an extra accuracy. „_., „„„ 

Both of these problems have been overcome in |7Z'] (see also p3] ) for the case of the Navier- 
Stokes problem in a strip. The key idea there was to use the special weights 

ipeix) := (l + e^|x-xo|)"^ 

depending on a small parameter e and this small parameter is allowed to depend on the solution. 
Then, the extra cubic term ip'^u^ ~ sipu"^ can be made small by the proper choice of the parameter 
e depending on the "size" of the solution u. The second problem related with the pressure has 
been overcome by multiplying the equation by mp — Vip where v^p is a small corrector which 
makes this multiplier divergence free and which can be found as a solution of the proper linear 
conjugate equation, see p2j for the details. 

The main aim of the present paper is to extend the weighted energy method to the case 
of (damped) Navier-Stokes equations in the whole space. Note that the analytic structure of 
the pressure term is essentialb!^ different than in the case of a strip with Dirichlet boundary 
conditions, so the method of I{z2] does not work in a direct way and requires an essential modi- 
fication. In particular, we have to handle the pressure term in a completely different way based 
on the explicit formula for Vp through u®u via the convolution operator and on the trick with 
the proper splitting of the convolution kernel developed in \f^\ which allows us to treat Vp in 
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weighted and uniformly local spaces. Thus, combining the ideas from p^J and I'TB'J, we will give 

a comprehensive study of the damped Navier-Stokes equations in the uniformly local spaces. 

The main result of thejjaper is the global well-posedness and dissipativity ofllie damped 
Navier-Stokes equations (pi.ip in the uniformly local space L^(M^) (see Section |E] tor the def- 
initions). Namely, if g £ L'^{M?) is such that divg' = and curler G L°°(R^), then for every 
divergence free uq € LI(M'^) there exists a unique global weak solution u{t) G L^iM?) of (pi. Ill and 
the following estimate holds: 

'^' + Qi\\9\\a + \\cnTlg\\L^), t>0. 



(1.5) 



u 



<Q{\\uo\ 



Ll)e 



where positive constant /3 and monotone increasing function Q are independent of t and u. This 
dissipative estimate allows us to apply the highly developed general theory of dissipative PDEs 
for the further study of the damped Navier-Stokes equation (jll.ip . In the present paper, we 
restrict ourselves only to verifying the existence of the so-called locally compact global attractor 
for this problem, further applications will be considered in the forthcoming paper. 

Although the most part of the paper is devoted to the damped case a > 0, we also consider 
the applications of the developed theory to the classical Navier-Stokes problem where a = 0. In 
this case, as elementary examples show, the solution u{t) may grow at least linearly as time goes 
to infinity, so the long standing open problem here is to obtain the realistic/sharp upper bounds 
for the growth as f — )• oo. Although we did not solve this problem conipletely, the developed 
method allows us to improve drastically the best known upper bounds (pi.4|) and to show that 
the solution u{t) can grow not faster than quintic polynomial in time: 

(1.6) ||n(t)||^2<C(t + l)^ 



where the constant C depends on uq and g, but is independent of t. Moreover, we obtain optimal 
bounds for the growth rate in time of the mean values of the solution u{t) with respect to the 
balls of time-dependent radiuses R{t) := (t + 1)^: 



\u{t,x)\'^ dx <C{t + l), xqG 



(this quantity indee 
tions), see Section] 



R{ty^ J\x-xo\<R{t) 

i grows linearly in time for simple examples of spatially homogeneous solu- 

...^,, ^^^ ^^^^.^^. i^foT more details. , . ^ 

The paper is organized as follows. In Section we recall the definitions and basic properties 
of the weighted and uniformly local Sobolev spaces, introduce special classes of weights and 
derive a number of elementary inequalities which wUl be used throughout of the paper. 

In Section 1^ we introduce (inspired mainly by l [iaj ) a number of technical tools which allows 
us to treat the pressure term in the proper weighted and uniformly local spaces as well as to 
exclude it from the various weighted energy estimates. 

Section H^^is devoted to the case of infinite-energy solutions which however decay to zero as 
|x| — )• oo (although the rate of decay can be arbitrarily slow). In this case, we are able to verify 
the global in times boundedness of solutions based solely on the weighted energy estimates 
(without using,the vorticity equation (|li.2|) ). 

In Section IQ" we treat the general case of spatially-non-decaying solutions u{t) G -L^(]R^). The 
global in time bounds for such solutions are obtained combining the weighted energy estimates 
obtained in Section IfflVith the vorticity estimate ()li..3|) . . 

In Section |Q^ we verify the dissipative estimate (|ll.b|) and show that the solution semigroup 
associated with equation (jll.l|) possesses a locally compact global attractor. 

In Section 0^ we consider the classical Navier-Stokes problem (a = 0) and derive thejsolyno- 
rnial unper bounds for the growth in time of the solution u{t). In particular, estimates (jll.6|) and 
(jll.V|l are verified here. 

Finally, some approximations of divergence free vector which are necessary in order to prove 
the existence of solutions in uniformly local spaces are considered in Appendix 1 (see Section Qj' 



and the proof of one important for our purposes interpolation inequality is given in Appendix 2 
(see Section r^ 
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2. Preliminaries I: Weighted and uniformly local spaces 

In this section, we briefly discuss the definitions and basic properties of the weighted and 
uniformly local Sobolev spaces (see I'M f i^Z f ilMftor more detailed exposition). We start with the 
class of admissible weight functions and associated weighted spaces. 

Definition 2.1. A positive function 0(x), x £ M?, is a weight function of exponential growth 
rate // > if 

(2.1) ^{x + y)<Ce^'^y^(p{x), x,yeR^. 

The associated weighted Lebesgue space L^(M^), 1 < p < oo, is defined as a subspace of functions 
belonging to L^^^(M^) for which the following norm is finite: 



IP 



{x)\u{x)\^ dx < oo 



rl,P 



and the Sobolev space W, (M ) is the subspace of distributions u G P'(M ) whose derivatives up 
to order / inclusively belong to L^(M^) (this works for positive integer / only, for fractional and 



rl>P 
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negative I, the space W^ is defined using the interpolation and duality arguments, see \B, 22] 
for more details). 



The typical examples of weight functions of exponential growth rate are 



-e\x—xo\ 



or 



-y/l + £^\x-Xo\^ 



e G 



xo G 



Another class of admissible weights of exponential growth rate are the so-called polynomial 
weights and, in particular, the weight function 



(2.4) 



'XQ 



(x):= 



1 



1 + Ix — Xol 



xo G 



which will be essentially used throughout of the paper. 

Next, we define the so-called uniformly local Sobolev spaces. 

Definition 2.2. The space L^{E?) is defined as the subspace of functions of L^^^(M^) for which 
the following norm is finite: 

(2.5) 



\U\\tP 



sup ||n||iP(Bi ) < oo 

X0&? ° 



(here and below B^ denotes the i?-ball in M^ centered at xq) and the space L^(R^) is a closed 
subspace of L'^{M?) which consists of functions tending to zero as |x| — )• oo in the following sense: 

(2.6) ^ lim \\u\\lp,bi\ =0. 



lim ||n||ip(m ) 

|a;o|^-co 



rl,V( 



l^) (resp. W^'^{E?)) are defined as subspaces of distributions u G P'(M^ 
whose derivatives up to order / inclusively belong to L^(M-^) (resp. L^{M?)). 



The spaces W^ 

Pm2^ ^.ggj. rPm2^ 

In slight abuse of notations, we also define the spaces Wj'^([0,T], 1^™'''?(R^)) of functions 
depending on space x and time t variables using the norm: 

(2.7) 



sup sup \\u\\Y^l,pl\Qrp-,y^rn,qlBl\\<00. 

te[0,T-l]xoGR2 "° 



IPIIV1/''P{[0,T],VK™'9) '■- 

The next proposition gives the useful equivalent norms in the weighted Sobolev spaces 



Propl . equiv 



Proposition 2.3. Let (p be the weight function of exponential growth rate and let 1 < p < oo, 
/ € R and R > 0. Then, 



1 . equiv 



(2-8) Ci / (t>{xo)\\u\\^i^p,j,dxo<\\ur,^j,<C2 (/>(xo)||n||^,,j,.^ dxo, 

where the constants Ci depend on R, I and p and the constants C and fi from iU.iW . but are 
independent of u and of the concrete choice of the weight (j). 



.eight 
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For the proof of these estimates, see e.g., I[5]. 
Thus, the norms f 



2^0 fc 



[xo)\\u\\\vi 



, p.„jj . dxo computed with different iZ's are equivalent. 



Propl .more 



The next Proposition gives relations between the weighted and uniformly local norms. 

Proposition 2.4. Let cj) be the weight of exponential growth rate such that J^ ^.2 4>dx < oo. 
Then, for every u £ W,'^{M?) and every k > 1, 



1 .ul-w 



(2.9) 



\u\ 



W''P{B-) 



< C 



y<^B- 



^\\w'-p(Bl) dy<C^ ^ 4>{y - 2;o)lkll^!,p(Bi) dy 



j/eiR2 



l.ul-wl 



and, in particular, fixing k = 1 in (||2.9|I and taking the supremum with respect to Xq E 
have 

(2.10) 



we 



■u||r^i,p < C sup ||n| 



4i{--Xq) 



1 . w-ul 



where C is independent of u and the concrete choice of the weight (p. In addition, 

(2.11) Ml,,<c\mL4u\\l,,, 



w:: 



'wt 



where C is also independent of u and the concrete choice of (p. 

1707^703 

For the proof of these results, see e.g., \p2^ i24j . 

The next lemma gives a simple, but important estimate for the weights Oxoi^) which will 
allow us to handle the convolution operators in weighted spaces. 

eta 



Leml. weight Lemma 2.5. Let 9xo{x) be the weight defined via (jl2.4l) . Then, the following estimate holds: 

l.thetakeyl (2.12) / 6*^,0 (a;)6'j^o (x) dx < C6l^o(yo), 

' JxgR2 

where C is independent of XQ,yQ £ M?. 

Proof. Indeed, due to the triangle inequality, 

1 1 

2 + \x- xo\^ + \x- yo\^ > -(1 + (|x - xol + \x - yo\f) > -(1 + |xo - yo|^)- 

Therefore, 

^.•o(a;)^j,o(^) = 2 + |x-xo|3 + |x-yoP Vl + |x-xo|3 + l + |x-yolV " 

and the integration of this inequality with respect to x gives the desired estimate 
the lemma is proved. 

Corollary 2.6. Let Oxoi^) be defined via (||2.4|) . Then, for every u £ L^ C^"^), we have 




Corl . conv 



1 . inv 



(2.13) 






(pyo{xo)\\u\fp dxQ, 



1 . convol 



1 . convest 
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, R-thetakey 



where C is independent of hq € M. Moreover, if the function v G ifoc( 
(2.14) \v{y)\<C f ey{x)\u{x)\dx, Vy G M^ 



satisfies the estimate 



then V e L^ (M^) and 



(2.15) 



If II rP < Cl ||n|| rP , 

"vo vo 



where Ci depends on C , hut is independent of yo £ 



Pj'ooJ. Indeed, ()l2.i3|l is an immediate corollary of ()b.l2|) and the Fubini theorem. To verify 
(jl^.ib ). we note that, due to the Holder inequality, estimate (jl^.Mp implies that 



HyW<C ey{x)\u{x)\'P dx. 

JxeR'2 

Multiplying this inequality by Oyg{y), integrating over y and using the Fubini theorem and 
we have 



takey 



IP 



iwir.p < c 



"^0 Ja;GR2 Jy^M.'^ 

and the corollary is proved. 



yoiy)OAy)dy\uix)\Pdx < Cl I ey,{x)\u{x)\P dx = Ci\\u\\l, 



D 



We conclude this section by introducing some weights and norms depending on a big parameter 
R which will be crucial for what follows. First, we introduce the following equivalent norm in 
the space W/'^(M2): 

(2.16) 



I^IIh/^.p 



sup \\u\\wi,p/BR)- 



Then, according to (|l2Tn|? 

(2.17) \\u\\^i,, < \\u\\^i,, < CR^/P\\u\\„,i,: 

where the constant C is independent of i? > 1. We also introduce the scaled weight function 



iiy^.P' 



(2.18) 



i,XO 



[x):- 



1 



R^ + \x — xo\ 
Then, the scaled analogue of (|l2.ij|l reads 



R-'e,g/R{x/R). 



(2.19) \\nr^,,.^^n) < CR-' 



W'.P(B«) 



dy < C^R 



^R,x„{y)\\u\fy^^,^p^j^R^dy, 



[1 t.lie 
where the constants C and C^ are independent of R and the scaled analogue of (Il2.12|l 

(2.20) / OR,,,{x)eR,y,ix)dx<CR-^9R^,,{yo), 



takey 



where C is independent of xo,yo S 
dR,yo{xo), integrating over xq and using 



and„i?,> 



(ILl.^ul). we se( 



11 iil-uR 

0, Moreover, multiplying inequality (jb.lujl by 



see that 



l.eqtheta (2.21) 



xeR^ 



:,xo 






dx. 



where C^ is independent of R. We also note that, analogously to (|l2.1l|) and using ifz.Li 



ulR 



(2.22) 



xeR2 



7R,XQ 



\ Ml llw''P{BR) - II llw''P - ^ II llw 



where the constants C and Ci are independent of i? ^ 1. Finally, we need one more simple state- 

h]1 R • ] ^ n 



ment about the behavior of the quantities involved into (jl2.22p in the case when u G Wj^ (M ) 



rl,Pfm,2\ 



1,P{ 



Propl. decay Proposition 2.7. Lei u G T4^^'^(M"'). Then, 



limzero 



gdiv 



0.5 



0.6 



. conv 
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0.7 



0.8 



0. Igood 



1 



(2.23) i5Li?2ll^ll<-(M^) = 0- 

The proof of this proposition is straightforward, so we leave it to the reader. 

3. Preliminaries II: Excluding the pressure 

In this section, we introduce the key estimates which allow us to work with the pressure term 
Vp in the uniformly local spaces. Note that the Helmholtz decomposition does not work for the 
general vector fields belonging to L^(M^), so the standard (for the bounded domains) approach 
does not work at least directly ^nd we need to proceed in a bit more accurate way. 

As usual, we assume that (jILlJ) is satisfied in the sense of distributions. Then, taking the 
divergence from both sides of ()ll.l|) and assuming that the external forces g are divergence free: 

(3.1) divg = 0, 
we have 

2 

(3.2) A^p = dw{{u,Vx)u) = ^ dx^dxjiuiUj). 

Thus, formally, p can be expressed through u by the following singular integral operator: 

(3.3) p(y) ■= '^Kij{x-y)ui{x)uj{x)dx, Kij{x) ■- '^ 



^tL"itIy ■j 



27r 



which we present in the form 

(3.4) p = Kw := K *w, 



w := u (S> u, K * w = y^ Ki 



*Wi 



«j 



It is well-known that the convolution operator K is well-defined as a bounded linear operator 
from w G [L'^(IR^)]'' to p G L'^(Mp), 1 < q < cx), but it is not true neither for q = oo nor for the 
uniformly-local space L^(R^). However, as the following simple lemma shows, the gradient oi p 
(which is sufficient in order to define a solution of (lli.i|l ) is well-defined in uniformly- local spaces 
and has natural regularity properties. 



Lemma 3.1. The operator w — )■ VxP, where p is defined via (113. 4|) can be extended by continuity 
(in Lq^J in a unique way to the bounded operator from [L^(M^)]^ to \W^ ''^(M'^)]^, 1 < g < oo. 

Proof. Indeed, let ipix) G C^(M^) be the cut-off function being identically one near x = (for 
|x| < 1/4) and identically zero for |x| > 1/2 and let w G [C^{^^)Y. Then, 



(3.5) Vyp = yyY, I Hx- y)K,j 



(x — y)wij{x) dx+ 



+ y^ V„[(l - ip{x - y))Kij{x - y)]wij{x) dx := Kiw + K2W. 



The first operator is extendable to the uniformly local spaces since^the kernel vanishes if |a;| > 1/2 
and by the standard Zygmund-Calderon estimates (see, e.g., ptTjjt we have 



(3.6) 



where C is independent of.xn G M^. Thus, multiplying (||3 61) by d 
and using (||2.9|) and (||2.»|) . we see that 

(3.7) ||lCiti;||„,-i,, < C 



[xq), integrating over xq G 



Wl 



\w\ 



U' 



0. Iggood 



0.2ggood (3.9) 



0.9 
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holds for every weight function of exponential growth rate. Moreover, taking (j) = (J)xq{x) := 
' (x—xq) in (jl3.V|l where ia ari integrahle weight of exponential growth rate, taking the supremum 



"S 



over xq G M and using (|l2.iup and (p. lip , we see that the operator Ki can be extended in a 



unique way on functions w G [L^(M^)]^ with the natural estimate 
(3.8) ||Kii(;||„,-i,9 < C||w||r9 

\ J II J- IIW — II lli^fc 



The second operator of 



is a smoothing convolution operator with integrable kernel 



(decaying as ^3- as |3;| — )• 00), namely. 



\K2{z)\ < 



C 



1 + |Z|3 

for some constant C independent of z. Therefore, according to Corollary 



,1 . conv 



'xo 



(x):- 



where C is independent of xq £ 1 

can also be extended by continuity to u; G [L^(M^)]^ and 



and the lemma is proved. 



l + |x-xo| 

Thus, taking the supremum over xq G 
li 



^, we see that K2 
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We will denote the operator obtained in the lemma by V^-P and the corresponding term in 
the Navier-Stokes equation will be denoted by VxPiu u). Then, in particular 

(3.10) 



I)2m2 



VxP{u®u):[Ll{^')] 



\W'. 



-1.2/ 



e)]2. 



We are now ready to give the definition of a weak solution of problem 

^) be divergent free. A function u{t,x) is a 



Definition 3.2. Let the external forces g £ Lf^ 
weak solution of problem (jll.lj) if 

(3.11) nGL~([0,T],L2(M2)), V^u € lI{[0,T] xR^) 

and the equation is satisfied in the sense of distributions with VxP = V^P(u (X" u) defined in 
Lemma II3. II Note that, according to the embedding theorem, u G Lf{[0,T] x M^) and u u £ 

a3,r] X M2) and, due to the previous lemma, V^P G Ll{[0,T],W^'^'^(M?)) and the equation 
can be understood as equality in this space. 

Remark 3.3. We emphasize once more that only the gradient of pressure Vp is well-defined as 
an element of L^([0,T], W^ ' ), but the pressure itself may be unbounded as |x| — )• cx). To be 
more precise, the operator VP defined above satisfies 



div VP(w) = ^ dxi dx^ Wij , curl VP{w 







in the sense of distributions. These relations can be justified by approximating w by finite 
functions and passing to the limit analogously to Lemma ll3.ll Therefore, there is a function 
p G L'^{[0,T],Ll^{R^)) such that 

Vp = VP{w), 

see Ifib], but this function may grow as |x| — )• 00. ISlpte also that the choice of Vp = VP(n u) 
is not unique. However, if pi and p2 both satisfy ()l3.2p (for the same velocity field u), then the 
difference pi — p2 solves A(pi — P2) = (in the sense of distributions) and, consequently is a 
harmonic function. Moreover, every harmonic function with bounded gradient is linear, so, if 
we want the velocity field u to be in the proper uniformly local space, the most general choice 
of the pressure is 



(3.13) 



Vp = \IP{u® u) + C{t), 



Rem.non-div 



g-proj 



0.11 



0.12 
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where the constant vector C{t) depends only on time (and is independent of x) and VP is 
defined in Lemma IkTi In the present paper, we consider only the choice C{t) = 0. In a fact, 
the vector C{t) should be treated as one more external data and can be chosen arbitrarily, but 
this does not lead to more general theory since everything can be reduced to the case of C = 
by replacing the external force g hy g — C{t). 

bd i y 

Remark 3.4. Assumption (|ia.ip on the external forces g can be slightly relaxed. Namely, we 
may assume instead that 

(3.14) 

where 11 is the prop 
this case, instead of 



Ug G [LUr')] 



2m2 



tension of the Leray projector to the divergent free vector fields. In 
, we have 



AxP = ^da;,d. 



■j {^{"^j J 



divg 



J, J 



which gives 

V^p = V^P{u ^u) + V^.(A^.)-i divg = V^P(n ® n) - (1 - n)^ 

1 0.1 

and, inserting this expression to (|li.l|) . we see that the problem is factually reduced to the case of 
divergent free external forces by replacing g by Ilg. However, even the assumption g G [L°°(]R^)]^ 
(or [C7~(m2)]2^ 



as log |x|, see 



does not imply that Ug £ [L'^(M^^)]'^ . In this case, we can only guarantee that 



but the functions with bounded mean oscillation may grow as |x| — )• oo (say. 



]), so 



IS not satisfied automatically. 



{u, Vx)u = AajU — au + VxP{u <S> u) + Ug 



Note also that this condition is necessary if we want the solution u{t) to be bounded as 
|x| — )• oo (and belong to the proper uniformly local Soboloev space). Indeed, after replacing the 
pressure term, (jll.i|) reads 

dtu 

and this equation cannot have bounded (as |x| — )• oo) solutions u{t) if Ug is not bounded (as 

l^l — )• oo). 

We conclude this preliminary section by establishing the key estimate which allows us handle 
the pressure term in weighted energy estimates. To this end, we introduce for every xq G M^ 
and i? > 1 the cut-off function ^r^xq which satisfies 



fR,xo{x) = 1, for X e B^ 'PR,xo{x) = 0, for x ^ B.^ 



2R 

Xo 1 



(3.15) 

and 

(3.16) \Vxm,x,{x)\ < CR-'ip](l^{x), 

where C is independent of R (obviously such family of cut-off functions exist). Then, the 
following result holds. 

Lemma 3.5. Let the exponents 1 <p,q <oo, - + - = 1, w G [L^(R2)]4 ^nd v G [^^^'^(M^)]^ be 
divergence free. Then the following estimate holds: 

(3-17) \('^xP{w),^R,xoV)\<C / 0R,xo{x)\\w\\Lp(^BR)dx-\\ip](lv\\Li, 

where C is independent of R and xq and Or^xq (x) is defined by ()l2.18|l . 

Proof. Indeed, it is sufficient to verify the formula for w,v £ C^ only. For general tf's and u's 
belonging to the uniformly local spaces, the estimate can be justified in a standard way using 
the approximations. For the functions with finite support, the function p is well-defined, so we 
may integrate by parts and write 

(3.18) C^xP{w),ipR,xo1') = -{P{w) - Pxo,'^x^R,xo ■ v), 



10 

I P 7 

where px^ = Pxoiu, R) is a constant which wih be specified below. Analogously to (|l3.b|) . we write 
0-l6| (3.19) p{y) = V / i^nix - y)Kij{x - y)wij{x) dx+ 

+ Xl / [(1 ~ ^^(^ ~ y))Kij{x - y)]wij{x) dx := Ki^r{w) + K2^r{w). 

where iPr{x) = ip{x/R) and ijj is the same as in {p.b\) . Actually, (|b.lijp coincides with lp.b\] up to 
scaling x and absence of the gradient. Thus, using the scaling x — >• Rx, xq — )• Rxq and y ^-t- Ry, 
we obtain the operators independent of R. Estimating the first one analogously to (IH.b|) and 
returning back to the unsealed variables, we have 

(3.20) ||^i,ij(ty)||LP(ij2fl) < C||u;||ip(B3fl) < 

<CiR^'^ \\w\\LPmii)dx <C2R 9R^xo{x)\\w\\Lp,BR)dx, 

J\x-xo\<AR "" JR^ 

where we have also implicitly used fc.iu|l . 

The second operator K2 r is a bit more delicate. Indeed, in contrast to (jl3.i3|) . the convolution 



0.17 



0.18 



0.19 



0.20 



kernel is not differentiated by y and decays only as -r^^ as |x| — t- 00 which is insufficient for 



the uniformly local and weighted estimates. However, following p3], we use the trick with the 
properly chosen constant p^^ which allows to resolve the problem. Namely, we fix 

(3.21) pxo ■■= K2.r{w){xq) = y^l (1 - ^Rix - xo))Kij{x - xo)wij{x) dx. 



Then, the new convolution operator 



K2,r{w) -Pxo = ^ Kij{x,y) 



Wij dx 



has the kernel 

Kij{x,y) := (1 - iI)r{x - y))Kij{x - y) - (1 - iJjr{x - xo))Kij{x - xq) 

which decays as \x — y\~'^ when |x — y| — t- 00 and |y — xqI remains bounded, so the singularity at 
infinity disappear and the uniformly local estimates work. Indeed, in scaled variables x = Rx' , 
y = Ry', Xq = Rx'q, we obviously have 

(3.22) \Kij{x',y')\<C—j^ ^, \y'-x'o\<2 

1 + \x' — y'l'^ 

and, therefore, again in scaled variables 

\\K2,rH - PxohpiBl, )<C T-TZr—JTsMLr'iBl,) dx' 

and returning back to the non-scaled variables, we see that 

(3.23) \\K2^r{w)-Pxo\\lp{bIR)<C2R I ^3-n \^\\M\Lp{B§)dx. 

Combining (jOUj) and (I03J) . we see that 



w) -Pxo\\lv{bIR) < CR / GR,xo{x)\\w\\Lp(^BR)dx 

3f R and xq. Fina|lly. applying the Holder inequali 
SI) , we arrive at the desired estimate (|b.lY|) and finish the proof of the lemma. D 



11 . JTi tpart 

;e C is independent of R and xq. Fina|lly^pplying the Holder inequality to (|b.l?5|) and using 



11 



s3 



4. Global well-posedness: the case of spatially decaying initial data 

In this section, we will study the damped Navier-Stokes equations under the additional as- 
sumptions that 



3.gu| (4.1) 5GL^(R^), noGLj 

The main result of this section is the following theorem 



')■ 



ThO.l 



3. global (4.2) 



0.13 



. huuge 



2.13 



Theorem 4.1. Let the external forces g and the^initial data uq be divergence free and satisfy 
()Kl.i|l . Then, the damped Kauicr- Stokes problem (jli.i|) possesses at least one global weak solution 
in the sense of Definition J^J^ and this solution satisfies the estimate 



HT)\U 



+ l|Va.n||i2([o^-r]xR2) < Qiuo,9), 



where the constant Q{uQ,g) is independent ofT > 0. 

Proof. We give below only the derivation of the key a priori estimate (|kl.2|) . The existence of 
a solution can be obtained after that in a standard way by apnroximating the infinite-energy 
solution by the finite-energy ones and passing to the limit, see I [i4ll'22j for the details (see also 
some more details at the end of the proof). 

We multiply equation (jli.i|) on ipr^xqU where R is a big positive nurnber which will be fixed 
below, xq E M^, and the cut-off functions ipR,xo satisfy (jp.lbp and (|k Ibp . Then, after the 
integration over x and standard transformations, we get 

(4.3) ~\Ht)\\h +a\\n\\h + W"^ A^Zo^)\\h ^ 



2dt 



VR 



Vfl.lQ 



<C\\gf^^ +C7i?-2||n||2 +Ci?^i||n||33 2^ +|(V^.P(n),(^^,^„n)|. 

Here the constant C is independent of R. To estiniate the term containing pressure, we use 
Lemma Ib.bl with g = 3 and p = 3/2. Then, due to (p.lYp together with the Holder and Young 
inequalities, 

(4.4) \{VxP{u),<fR,xoU)\ <C 0Rxo{x)\\u^u\\^3/2(BR)dx-\\u\\L3,B^R.) 



< 



<c 



<C{ I eR,xoix)dx 



lxoix)\\u\ 



LKBR) 



dx ■ \\u\ 



LHsm 



< 



1/3 



• 
\ 2/3 



OR,xo{x)\\u\\l^BR^dx\ • ||n||^3(B2H) 



< 



< CR~^/'^ 



2/3 



')r,xo{x)\\u\\^3^bR) dx] • ||n||i3(B2H) < 

< C* / 0R^XQix)\\u\\l3(^BR) dx + Ci?"^||u|||3(B2fl), 



where all constants are independent of i? ^ 1. Thus, 
d 



now reads 



(4.5) 



0.25 



dt ""v'-R 

<C\\g 

We introduce 

(4.6) 



U{t)\\l2 +q||u||^2 +\\VxU\\io.r.R.< 



^R,x 



\LHB^^) 



2 +CR \\u\\j2(f>2R\+CR I 



^lli3(B|fl)+'^ 



'^R,xo{x)\\u\\^3^^BR) 



dx. 



Zr 



S/O 



(n) := / 



:,i/o(a^o)hlli2 



dxn 



•PR^XQ 



12 



1 1 .fiqt.heta 

Then, using (jl2.2i|) . we have 



3.Z 



(4.7) 



Co 



■,yo 



u 



\l2,BR)dy< ZB^,yo{u) <Ci 



■,yo 



u\ 



L^BR) 



dy, 



where Cj are independervt otR- Multiplying now equation fi.b\] on 9fi^yg{xo), integrating over 
xo £ M^ and using (jl2.2U|l . we see that, for sufficiently large R, 



3. main 



(4.8) -Zfi,,,,(n(t)) + 2f3Zn,.Mt)) + 2/3 



:,a;o(^)lFllvyi.2(_BH) dx < 



< CZR^^,{g) + CR ^ Or.xo{x)\\u\\Iz,j^r. dx, 



3. int 



where the positive constants C and /3 are independent of R. 

Thus, we only need to estimate the integral containing the L'^-norm of u in the right-hand 
side. To this end, we use the interpolation inequality 

(4.9) Il^lli3(ijfl ) < C'll^lli2(iJfl )hllwl>2(BB ), 

where C is independent of xq and R. Then, using also the Holder inequality, we end up with 



3.qubic| (4.10) CR ^ eR^^^{x)\\u\\ls,^n.dx<CR ^ OR^^^{x)\\u\\l^,^n)\ 

' JxeR'2 ^ '"' Ja:eR2 ^ "^ ' 



U U/1,2 



W^'HBR) 



dx < 



< CR-'^Wul 



f 2 



■,xo[X)\m\L2{Bri)\m\w^''^iBR) 



dx < 



<CR'^\\uur2 



Zr,xo{u) +13 9r^xo^ 



X U 



W/i.2(Bfl) 



dx. 



0.27 



Inserting this estimate into the right-hand side of ()I4.8|I . we arrive at 
(4.11) ^ZR,,Ju{t)) + pZR^^^{u{t)) + /3Zfi,,„(n(t))(l - CR-^\\u{t)y^ 

at b.R 



U )+ 



+ /3 / 6'ij,^-o(a:^)lkllvi/i>2(Bij) dx < CZR^^^ig)), 



where the positive coristants C and /3 are independent of R. 
Finally, due to (jy'iyll " 



0.281 1 (4.12) 



-21, 



3. last 



\\u\\l2 < CR sup ZR^^,g{u), 
and (|H?TT|) reads 
(4.13) -^Zfi,,o(n(i)) + /3Zij,,o(n(t)) + /3Z/j,,„(n(t))(l - KR-' sup Zfi,,,(u(t)))+ 

+ 13 ^R,xoix)\\u\\'^i,2,BR) dx < CZR^^f,{g)), 

where the positive constants K, (3 and C are independent of i? ^ 1. 



e claim that the key estimate (|kl.l3|l is sufficient to finish the proof of the theorem. Indeed, 



.13|) gives the desired estimate 

rt+l 
0T28] (4.14) ZR,,„{u{t))+ II 



dt < 



< CZR^,,iuo)e-^' + CZR,^,{g) < C{Zr,,,{uo) + ZR^^,{g)) 
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0.29 



3.un 



Th3 .unique 



uniqueness 



2.unieq (4.20) 



obal 



(which imphes (jl4.2|l ]ust by taking the supremum via xq G 
in such way that 



if we are able to fix i? = R{uq, g) 



(4.15) 



R- 



Ht)) < 



1 

K 



sup Z/j,^o 

xqGK 

holds for every t G M^. 

On the other hand, due to the decay assumption uo,g G L 
fix i? = R{u(),g) in such way that 



-——. 1 .decay 
and Proposition Il2. '4 we mav 



0.30 


(4.16) 
where C i 

and p.UI 

(4.17) 






) hoi 


3 . cont 





CR-\1 + Zr,^,{uo) + ^i?,xo(5)) < 
K are the same as in 



1 
2K' 



xo G 



Hj) and (Ikt.l^jj) respectively. ^^^ 

is satisfied for t = and the standard continuity arguments show that both (JKJ.ibp 
hold for all t G M+ if we a priori know that the function 



t 



sup Z/J,xo 

a;oGK2 



(n(t)) 




obal 

IS verified under the 



is continuous as a function of time. Thus, the 
additional continuity assumption for the function 

However, in contrast to the case of finite energy solutions, the regularity of a weak solution 
stated in Definition Ikj"^ I is not enough to establish that u € C([0, T], L'^{M?)), so a priori, function 
(jK-iy]) may be not continuous (in a fact, it is continuous a posteriori due to the decay assumption 
uo,g G L^(M^), but in order to verify that, we need some extra estimat^- 

We overcome the continuity problem in an alternative way (following \ [Z2\ ) , namely, we cut-off 
the external force g and consider the&fiie support external forces g^ (x) := g{x)(pNfl{x) (where 
'Pr,xo{x) are defined by fe.lbp and (JD.16P ). approximate the divergence free initial data, wn b 



the sequence Uq of the divergence free vector fiejd^ with finite supports using Corollarv lfe.^i see 
Appendix, and will apply estimates (||4.14|) and (|l4.1bp only for the approximative solutions u^{t) 



^ 



N 



u^ € C([0,r],L^(M^))) for which the continui 
31 we may fix i? = R(uo,g) ^ 1 in such way that 




,N 



which correspond to the initial data Uq and external forces g . Indeed, for every A^, u is a 



finite-energy solU||40 
due to Corollary 



IS obvious and, 



will hold uniformly 



for all N. Thus, we have justified the uniform with respect to N estimate 

r-t+l 

(4.18) ZH,,,(n^^(t)) + 



.A^ 



Wi 



dt < C(Z/j,^o(no) + Zij,^o(c/)), 



p2 



where C is independent of i G M+ and xq G M'^. Passing now to the limit A — >■ oo (in the local 
topology of -L^^^(M^J), we construct the desired weak solution u of the Navier-Stokes equations 
which will satisfy wi.2\) and finish the proof of the theorem. ~ 



D 



Remind that we have construe 
which satisfies the key estimate 



a weak solution u{t) of the Navier-Stokes problem 
Tip to the moment, we do not know whether or not all 



weak solutions satisfy it.^ Jlowever, the next theorem shows that the weak solution is unique 
and, by this reason, fL.'Zh holds for all solutions. 

I tp,e sense of Definition b.^ is unique. Moreover, for any 
we have 



Theorem 4.2. The weak solution t 
two solutions ui and U2 of problem 

(4.19) ||ni(t) 



U2 






IS unique. 

<Ct||^i(0)- 112(0)11^2 , 
where the constant C depends on T, R, ui and U2, but is independent of xq G M?. 

1 0.1 

Proof. Let ui{t) and U2{t) be two weak solutions of (|li.ip and v{t) := ui{t) — U2{t). Then, this 
function solves 



dtv + {ui,\7x)v + {v,\7x)u2 = AxV - av + \7P{wi - W2), divu = 0, 
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where Wi = ttj (8) jUj Multiplving this equation on ^r^xqV (where the cut-off functions ipR^xoix) are 
the same as in (|l3.1bp and (|l3.i6p ). after the integration by parts and straightforward estimates, 
we get 

lAlU,l|2_ ,.,IU,I|2^ , IIY7 /,„ 1/2 M|2 _ ^^||Y7 „Jk„._„„ IU„ 1/2 ||2 



3-^^^dif| (4.21) 1:^^11 +"II^IIl2 +||V^(?;(/j^/ )||^2(^2fl)<C7||V^n2||i2(B2fl)||wc/?^; 11^4 + 

+ Ci?"^||ni||i2(B2B)||f||^4(B2fl) + |(VP(u;i -W2),VtfR,xo)\ +Ci?"^||f||^2(B2fl). 

jr.fiTiip . 2 
We start with estimating the most complicated term containing pressure. Using Lemma 113. bl 

together with the fact that ||iij(i)||?2|'Rfl -i is bounded (by the definition of a weak solution), 
4]) , we have 



3. int-lad 



(4.22) \iVPiwi - W2),VipR^xo)\ < C / Or,xo{x)\\wi - W2\\L'i/'-ilBR)dx ■ \\v\\^i(^2R) < 

<C Or^xo{x){\\ui\\l2(^bR) + \\u2\\L^BR))\H\L'i{BR) dx ■ \\v\\l4(^b^^R) < 

J M 

< Cu^,U2 / SR,xoix)\\v\\Li(^BR) dx ■ \\v\\l4(^b^R) < 
Jr2 "- ^' ^0^ 

< c\H%iBi^) + cl„u,R'^ j^^ Or,x,{x)\\v\\1,^j,^^ dx < 

< C||u||^4(52fl) +CR OR^xoix^vWli^BR) dx, 

where the constant C depends on the L^-norms of the solutions ui and U2, but is independent 
of R. We now estimate the more standard first term in the right-hand side of (|hi.2i|) using the 
interpolation inequality 

(4.23) \\U\\%<C\\U\\l2\\VxU\\l2. 
This gives 

(4.24) C\\VxU2\\l2(bIR)\\v/rIJ% < Ci\\VxU2\\l2(b2^R)\HlI^^^^ INAv^rIo^Wl^ < 

|2 IU,||2 , ^lU, ^1/2 ||2 



< C2\\VxU2\\l2^b^^R-j\\v\\l2^^^^^ + l^WvfR^xoh^- 
tdif 

of Ui{t) are under the control, we end up with 



^'^R,xn 2 

Inserting the obtained estimates into the right-hand side of (|k[.2l|l and using that the L?-norms 



,)^ 



S.estgr (4.25) 3-||w||i2 - C\\VxU2\\l2(b2R)\\v\\l2 + W^ xv\\ r2 ( ^r ) 

< C'fllblli4(B2fl) + ^R I ^R,xoix)\Hl4(BR) dx, 



where the constants Cr depend on R and on the L^-norms of the solutions ui and U2, but are 
independent of xq G K^. Moreover, by the definition of a weak solutions. 



/ \\^xU2{t)\\l2(B2R) dt < Ct, 
Jo ^ ^0 ' 
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wh^^e^Qr depends on T, but is independent of xq. Thus, we can apply the Gronwall's inequahty 
which gives 



3. after 



(4.26) \\v 



2 



+ /3 / \\^xv{s)\\i2,B2R^ ds < CR,T\\vm\ 



2 



+ 



+ Cr^t \\v{s)\\l4(^B2R)ds + CR^T l9R^xoix)\\v{s)\\l4(^BR)dsdt 

where the constants /3 and Cr^t are independent of xn e ]R^._This allows us to multiply 
by ^_R,j;o(xo), integrate over xq G M? and use (jb.i^jl and (jl2.iu|l to obtain 



ter 



(4.27) \\v{t)\\l, +/3 / \\VMs)\\l2 ds< 



< Cr^t\\v{0)\\12 +Cr,t / OR^^„{x)\\v\\li,gn.dxds. 

A J a..aaj:-lad 

Estimating the L -norm in the right-hand side with the help of (|kl:.23|) . we finally arrive at 

/3 rt 



(4.28) \\v{t) 



lii 



+ 



R,^0 



\h 



\yxv{s)\\'h ds < CR,T\\vm\ 



2 



+ Cr. 



T 



\v{s)\\ 



2 



ds. 



JluLumeness 

Applying once more the Gronwall's inequality, we derive the desired estimate (j[4.i9|) and finish 
the proof of the theorem. D 

Remark 4.3. Note that the uniqueness is proved without assuming the spatial decay condition 
m.l\] and, therefore, holds for any weak solutions of the Navier-Stokes problem in the sense of 
Definition Ib.21 



s4 



Prop4. local 



4. local 



5. Global well-posedness: the case of spatially non-decaying initial data 

The aim of this section is to obtain the analogue of Theorem M.lJ without the extra spatial 
decay condition (jl4.1|l on the initial data. Recall, that estimate ()I4.2|) has been obtained using 
the purely energetic methods (weighted L^-estimates) without any use of the vorticity estimates 
(which are traditionally the key technical tools for studying the Navier-Stokes equations in the 
whole space M^). However, this energy method requires some decay of the initial data at infinity 
(although the rate of this decay may be arbitrarily slow) and we do not know how to remove 
this assumption remaining in the class of weighted energy estimates. 

In the present section, we will show that the extra decay assumption (|kI.I 1i £a|n be nevertheless 
removed if we combine the method presented in the proof of Theorem li4.ll with the classical 
vorticity estimates. Since these estimates require the initial data uq to be more regular than 
just uq G L^(R^), we start with reminding the results on ih^lqcal solvability and local smoothing 
property for the solutions for the Navier-Stokes equation 



mq, 



^\ii-^ ) he divergence free. Thei 
> and a unique weak solution of problem 



there exists time moment T = 
defined on the time interval 



Proposition 5.1. Let UQ,g E L^ 

nwuoWq 

te[o,T]. 

Proof. Indeed, the uniqueness of the weak solution is verified in Theorem 114.21 so we only need 
the local existence. To this end, as usual, it is sufficient to verify the proper local a priori estimate 
for the weak solution on a small timeinterval t £ [0, T]. To obtain such an estimate, in turn, it 
is sufficient to integrate estimate 
a supremum over xq S 

ft ft 

|2 I „ — / IIV7 „./'„Ml2 j„ ^ rilW^. I|2 I ||„||2 \ ^ ri \ IL./'„M|4 



i.llp (which holds for every weak solution) in time and take 
P. Then, after the straightforward estimates, we end up with 



(5.1) ||n(t)|| 



^b,R 



+ sup 

xoeK2 



|V^-u(s) 







liL 



ds < C{\\uo\ 



Ll 



b,R 



^b,R 



+ \\g\\j^ ) + C \\u{s)\\Uds 
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peal 



for some constant C which depends on R. Esthnate (|b.l|) is enough to conclude that there exists 
r = T(||no||i2 ,11511^2 ) such that 



b.R 



4. locest 



and the proposition is proved. 



D 



Prop4. smoo 



| Prop4. local 

Proposition 5.2. Let the assumptions of Proposition \b.l\ hold. Then, the local weak solution 
u{t) becomes smoother: u{t) G W^ ' (R^) for all t > and the following estimate holds: 



4. locsm 



(5.3) 



^L'^{[t,T],W^'^) 



l.2^ + \Mqi[t,T],w^-'') - ^^ ^^"^QiWaWq + IkoIlL^), 



|Pro p4. local 

where T <^ 1 is the same as in Proposition jit.ll t £ (U-^]- and the monotone increasing function 
Q is independent of the concrete choice of u{t) . 



Proof. To derive estimate (jb.3|L we multiply equation (|ll.l|) by t ^^ dx^{^R,xo9xiU) and integrate 
over X. Then, after the standard transformations, we get 



4. hi 



(5.4) 



d 
It^ 



\VxU\\\2 ) + 2t||A^u||^2 <C{t\\u\\L2^BlR)+t+'^)\W\\wi,2iB'^R) + 

+ Ct\iVxP{u<8) u,),Y,dx,{^R,xodx,u))\, 



where C may depend on R, but is independent of xq. Here we have implicitly used that in 
2D {AxU,{u,'Vx)u) = O^for every divergence free function and, therefore, the leading part of 
the inertial term in (jd.4|| disappears and the lower order remainder (which appears due to th£ 
presence of the weight) can be easily estimated by the first term in the right-hand side of ( lp.4|) . 
Thus, we only need to estimate the term containing pressure in the right-hand side of (jd.4|| . To 
this end, we observe that AxU is divergent free, so due to Lemma ll3.bf 



4. press I (5.5) Ct\{VP{w),ipR^xo'^xu)\ < Cit\\Axu\\L2 / 0R^xoix)\\u\\l4(BR)dx < 

1 "^-B.a^O JlR2 ^ ^ ' 

< t\\Axu\\l2^^ ^ +C2i leR,xo{x)\\u\\L2^BR)\Mw^.2(BR)dx) < 

<t||A^u||^2 +C3Q{\\uo\\Ll + \\g\\Ll) GR,xo{x)\\u\\lyi,2(BR) 



dx. 



sre the constants Q may depend on R, but are independent of xq G M^ (here we have also used 
) 2|) m order to control the L^-norm of u and the fact that t < T and Tis, small). Analogously, 
using that dxiU, i = 1,2 are both divergence free and using Lemma Id.5I with (pR^xo replaced by 
dx,^R,xo, we have 



4.pressl| (5.6) Ct\{VP{w),Vx'PR,xo ■ ^xu)\ < Cit\\Vxu\\L2(B'2R) Gr^xo{x)\\u\\14(bR) 



dx < 



< tlWrufr 



xu.\\L2(^B2R) ^^2 



+ C2 



i,xo{x)\\u\\L2(^BR)\\u\\w^.HBR) dx ] < 



< C||Vx.u||^2(B2fl) + CsQiWuoW^i + WgWq) / ^ Sr,xo{x)\\u\\'^i,2(^bR) dx- 
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4.hliiit 



4.glinf 



4. vor 



Prop4 . max 



4. max 



rop4. linf sm 



Th4 . main 



Inserting (jb.b|) and (jb.6|) into the right-hand side of equation 



and integrating in time over 



[0,t], we get 

(5.7) t\\VMt)\\h + I t\\^^ 
^«.-o Jo 



u[s) 



'^-«,.o 



ds < 



< C''3(II^oIIl2 + Mq) f / lh('S)ll^i,2(B2fl) 



+ 



^R,xo (x) 



<Qii\\uo\\Li + IIs'IIl2 



pcest 



where we have used (|p.2|) again in order to estimate the time integral of the VK^'^-norm of u. 
Since the last estimate js uniform with respect to xq £ M^, taking the supremum with respect 
to xq, we end up with lp.ii\] and finish the proof of the proposition. D 

We now make an extra assumption on the smoothness of the external forces g, namely, we 
assume that 

(5.8) curl<7 := d,,g2 - d,,gi G L°°(m2) 

and remind that the vorticity w := curlu satisfies the heat equation with the transport term: 

(5.9) dfUJ + {u,Vx)uj + au! — AxU = cniclg, a; | j^„ = curl uq • 

The crucial property is that the vorticity equation possesses the maximum/comparison principle 
which allows us to control the L°°-norm of the vorticity w. 

Proposition 5.3u..Let the assumptions of Proposition W.l\ tiolds~dnd let, in addition curluo G 
L°°(IR^) and Ip.iih is satisfied. Then, the following estimate holds for u :=curlu.' 



(5.10) 



\UJ 



-at 



1 



(pITf 



+ - curl 5 Lo 
a 



L°o < II curluollioo e 
is an immediate corollary of the comparison principle, see e.g. 



Indeed, estimate 
more details. 

Finally, if the initialxlata un G LI{M?) only as the next proposition shows, io{t) G L°' 
t > and Proposition lb. af can be nevertheless used. 

p4.1nra1 

Proposition 5.4. Let the assumptions of Proposition \\d^ ^l\ ^holds and let, in addition, 

'00/Tn,2 



rv] for 

e) for 



IHF be 



), oj{t) G L= 



for all t G (0, T] 



satisfied. Then, for any weak solution u{t) of problem 
and the following estimate holds: 

(5.11) ||a;(i)||ioo(]82) < t"^(5(||uo|lL2 + 11^11^2 + || curl^Hioo), 

for some positive N and monotone function Q which are independent of t, uq and g. 

Indeed, due to Proposition ISf t}(tj°°G Ll(M?) and u{t) G Wl''^(U?) for t > 0. The further 
regularity of w(t) can be now obtained by the classical smoothing estimates for the heat equation 
(p.9l) using, say, the Moser iterations. 

We are now ready to state and prove the main result of this section. 



Theorem 5.5. Letuo,g G L^(M^) be divergent free and let, in addition, curl (7 G L°^{E?). Then 
the unique weak solution u(t) of the Navier-Stokes problem (jll l|) exists globally in time and the 
following estimate holds: 



4.globalest (5.12 



\u 



Ll < Qdl^olli?) + QdblL? + II CUrl5r||Loo), 



where Q is independent of time and, therefore, every weak solution is globally bounded in time. 
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J Ei:Qp4.1inf sm 

Proof. We first note that, due to Proposition lb. 41 we may assume without loss of generality that 
curl no G L°°{M?), so we may use the vorticity estimate (JD.iup starting from t = 0. The idea of 
the proof is to estimate the L^-norm in the right-hand side of {fi.ii\) (see the proof of Theorem 
lii.ip m a better way using the vorticity estimate (p.iu|l and the proper interpolation. To this 
end, we need the following lemma. 

Lem4.key Lemma 5.6. Let the vector field u G [Wq' {B'^^)]'^ be such that divu, curl u G L°°{B'^^). Then, 

1/6 



4-keyest| (5.13) \\u\\l^bIR) < C||^II/2(bM) [\\ curl u\\ ^o. (^b^R) + II divu||icx,(B2H)J , 

where the constant C is independent of R and xq. Moreover, for any 2 < p < oo, 

1 1 



4.keyestl| {5 M) \\u\\l^^b^r-j < C\\u\\ ^2 ^b^r) {^\\ cml u\\LP(B2ii) + \\ div u\\LP(B'iR)j , ^ = 2 ~ 2(» - 1) ' 

where C may depend on p, but is independent of R and xq £M?. 

For the proof of the lemma see Appendix 2. ,„ , ^ 

We are now returning to the proof of the theorem. Using (|b.i3|) and the cut-off functions 
(jls.ibp and (jB.itip . we estimate the L^-norm of the solution as follows: 



, estcorrect 



(5.15) ||'u.|li3(5fl ) < \\uipr^xo\\l2(b^R) 



< 



5/2 1/9 

< C'll^llL2(B2i?) (II curl(v3/j,^.on)||L°= + || div {if r^xou)\\l^) ' < 
,5/2 /„„_!„ „ „„ „ \i/2 



<C\\u\\^,^^,^^[CR in||ioo(BM) + C||6j||i^(B2fl)j < 

^ rlU,l|5/2 II, ,l|l/2 ) rR-l/2||„||5/2 II ||l/2 

- '-^IFIlL2{B4fl)ll^llL=o -^ L'-rt ||W|li2(B4fl)|F||^oo(B2fl), 

where the constant C is independent of R and xq. Analogously, using estimate (|b.i4|) . say, with 
p = 4, we have 

1/3 / —1 \ 2/3 

II'"IIl°°(B2«) < \\V2R,xoU\\l°° < C'||n||^2(£4fl) [\\w\\L4(^BiR) +-R lkllL4(B4fl)j < 

< C\W\\w^'^{B*R) (-^ ll'^IU°° + R~^\W\\w^'^{B^R)) ^ 

< CR ' \\uj\\loo + C||n||^yl,2(54fl-), 

where the constant C is independent of R and xq and we have implicitly used the embedding 
theorem W ' C L . Inserting this estimate into the right-hand side of (jb.ibp . we arrive at 

f5 1fi^ Il7/I|3 <rll7/ll^/2 IL;ll^/2 +rR-l/2||„||5/2 II ||l/2 

)^<J.iUj \\"'\\l'^(BR ) — ^ ll"'lli;,2(^4ii-)||"^|lioo -|- on ||"'||^2/j54il-) ||"||^yl,2(g4ii') ^ 

^ r^U Il2 /ll ||1/2|| ||l/2 , n-l/2|| ii \ ^ 

^^\m\wh2(^BiR)[\ML°°\ML^BlR)^-^ W'^WL^BiR)) < 



lTyi>2(B«), 



< Ci?^/^ (i?"^||n||i2 + IIwIIl- ) ||n||^ 

11 .pq theta 

where the constant C is independent of R and xq- Therefore, using (|U.2ip . we have 



(5.17) / eR^x,,{x)\\u\\l3,j^R.dx<CR^/'^ (r '^\\u\\l2 +\\uj\\lo^) 6*^,^0 (x) ||n||^i,2(Bfl) ^^a;. 
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4. main 



4. loop 



4 . Zbound 



4 . bound 



4. best 



Inserting this estimate into the right-hand side of 
(5.18) j^ZR^,,,{u{t)) + /3Z^„,o(n(t))+ 



m 

^ITwe get 



+ (2/3-KR-^/^{R-^\ 



U\\r2 + \\U}\\fo 







R,xo 



[x)\\u\\l^,l,2^BR) dx < CZR^^^ig), 



where the positive constants C, K and /l^e independent of R and xq. 

Analogously to the proof of Theorem fe.il if the parameter R is chosen in such way that 

(5.19) KR-^/'^{R~^\\u\\l2 + ||w||ioo) < 2/3, 
the Gronwall estimate applied to (lb.is|i gives 

(5.20) ^i?,xo(0 < ^i?,.oK)e-^* + CZr^^M < CRiWuoWl 



+ \\9\\h] 



and, therefore, taking into the account (|b.lU|l and (jl2.19|l . we have 



(5.21) 



R 



-ii 



u 



^b,R 



+ lk(OllL°° < C'ldkolli? + II curluollioo + \\g\\^2 + II curl^llioo). 



Finally, arguing as in the end of the proof of Theorem Ikt.ll we show that both 
are satisfied if the parameter R is chosen in a such way that 



| 4-1nnn 1 4 .7h 

jbT9j and ()On 



Zbpund 



i?-V2 



(5.22) 

This, together with 
(5.23) 



/3 



uoll^,? + II curlno||L°° + IIs-IIl? + II curl^Hio 



KCi 

■ boimd 

.211) gives the estimate 



\u 



\u 



^b,R 



< C(||'Uo||^2 + II curluollLoo + ||gr||^2 + II curl 5-11^00 ) 

O 



and the theorem is proved. 



D 



s5 



6. DiSSIPATIVITY AND ATTRACTORS 

|S ^section, we have shown that the global weak solution of the damped Navier- 
exists and remains bounded as t — ?■ oo. The aim of the present section is 



In the previo 
Stokes system 

to derive the dissipative analogue of 
considered as the main result of the section. 



Do be more precise, the following theorem can be 



ThS.dis 



5 .dis 



5. wrong (6.2) 



Theorem 6.1. Let the assumptions of Theorem lb. d\ fiold. Then, the unique weak solution u{t) 
of the damped Navier-Stokes equations (jll l|) possesses the following estimate: 

-I3t 



(6.1) 



< 



o 



+ <5(IIS'IIl2 + ||curl5r||L° 



where the positive constant f3 and the monotone function Q are independent oft and uq. 
Rroqf. As in the proof 



Note also that 



em lb. b I we mav assume that curluo E L°°(M ) and use estimate 



.2U|| gives us the estimate of the form 

< CiWuo 



R~^\\ufr2 

^b,R 



' e-^' + \\g\\l2] 



^b 



which looks as dissipative. However, it is not sufficient to derive the desired estimate (Ib.lj) since 
the parameter R depends jin the initiaLdata, see ()D.22p . The idea of the proof of this theorem 
is to allow (following toJlUl^see also t^2l) the parameter R to depend on t {R = R{t)) and to 
derive the analogue of (|b.2|l for the time-dependent R satisfying proper dissipative estimate. To 
this end, we first need to know how the cut-off {(pR^xo) and weight (Or^xq) ^W'l^^'^s depend on 
the parameter R. We start with the cut-off functions. To satisfy (0.11) and (jp.lbp it is sufficient 
to take 
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where (/? is a single cut-off function which equals one at Bq and zero outside of Bq and satisfies 
(ll3.i6D with R = 1. Assuming that R = R(t) is a smooth function, the differentiation gives 



5.phiT (6.3) 



R{t) 

, <y / M ^ ^ P^ (^) r / Ml /2 



5 . thetaT 



S.maindis 



S.Rslow 



5. correct 



S.Rchoice 



where the constant C is independent of R and xq. Analogously, the straightforward calculations 
show that 

(6.4) 



|^Hw,.o(^)l<c^-W-^m^o(^)> 



where the constant C is also independent of R and xq. ,„ ■ 

Using these estimates and arguing exactly as in the derivation of (jb.i»|) . we end up with the 
following inequality: 

(6.5) ^^Zn^^^^^Mt)) + /3^ii(t),.o(^(i)) + (/3 - ^1 • ^^) ^i?,{t),xo(^)+ 



+ (2f3-KR{t)-^/\R{t)~^\ 



u\\t2 + ||w||roo) 



^R{t),xoix)\W\\^,^,^^Rit)^dx < 



<cz. 



,(9), 



'R{t),xo\y)^ 

where the positive constants C, K, Ki and /3 are independent of R and xq. 

Indeed, the estimates of terms which do not involve the time differentiation are identical 
and only the terms containing time derivatives may cause the difference. At the first step, we 
multiply (jli.ip by u^Pru\^^ and the term with time derivative now reads 



{dtu, uipR^t). 



XQ; 



1„ 



1 



and using 



iT 



"fRW.^o 



-{\u\ ,dtLpR(t), 



XOJ 



we estimate the extra term via 



\{\u\\dt^R^t),,J<C^ 



\R'{t)\ 

R{t) 



'L2(B-«W) 



One more extra term we obtain when we multiply equation (|kl:.b|l by 9R(t),xo ^^^^ integrate over 
Xq. Namely, 



d 



d 



XQi 



^i?W,.o(-o)^||«||,.^^^^^^^^xo = -Z«(,),,„(n) 



dtOR{t),yo{xo)\\u"^ 



Ll 



dxn 



'f'R(t),xo 



and the extra term here can be estimated using 

|2 



etaT 



Via 



^t^Rit),yo 



R' (f)\ 

(xo)||m||^2__ dxo\ < C p^^^ ZR{t),yQ{u) 



•fiR(t),xa 



R{t) 



The^above.two estimate leads in a straightforward way to the extra term in the left-hand side 

P-jLainais " •' 

of (jb.bj) . Ihus. under the extra assumption 



R{t) - Ki 

indis 



(6.6) 

(which we assume from now on to be true), (|b.fci|) will be identical to (|b.ib|) and we derive from 
it that 



(6.7) 

if R{t) is such that 

(6.8) 2j3-KR{t)-^''^{R{ty^ 



Rit)-^\\u\\r2 < C(\\U0\\r2e-^^/^ + ||5||r2) 



U 



L2 + W 

^b,R(t) 



L°° , 



> 
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5.RT 



5 . sem 



DefS.attr 



Cor5. attr 



5 . atrattr 



J 5.4Lorrect 

holds for all t. Note that, if (jl6.V|l holds, due to 



we have 



(6.9) R{t) 



-ii 



u 



L2 + W^ 



loo < C{\\uq\\i2 + II curlno||Ltx))e' 



-/3t/2 



+ IbllL? + llcurlc/llio 

D 

.1 



Thus, usinft^the continuitY arguments (exactly as in the proof of Theorem ll4.l|l . we may conclude 
that both (Ib.Vj) and (jl6.»|) hold if we take R{t) as follows: 



(6.10) 



R{t) 



|uo||l2 + ||curluo||L°°)e ''' + ||g||2,2 + || curl^Hioo 



where the parameter < 7 < (3/2 should be chosen small enough to satisfy 



\R'it)\ 

R{t) 



--ft I 



27- 



|no||i2 + ||curluo||L°o) 



|moIIl2 + II curlMo||L°°)e"T* + 11^11x2 + || curlfir||ioo 



1 

< 27 



low 

nhdeed, 



low 



and (llb.bll is satisfied for 7 = /3 mini i ^). 

Thus, estimate (|I6.Y|) is verified for R{t) satisfying (|b.iup and, therefore. 



B.fingal (6.11 



\u 



\u{t)\\ 



^b.R 



< 



< C'((lko|lL2 + ||curluo||L«=)e ''' + 11^11^2 + || curl 5(||loo 



FollLge 



-7* 



+ II5|Il2 



which gives the desired estimate (jl6.1|l and finishes the proof of the theorem. 
Thus, we have verified that the solution semigroup 
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(6.12) S{t)uo := u{t), S{t) : H^ -^ H,, H, := {u^ e [LI{^^)]\ divuo = 0}, 



where u{t) is ^ unique global weak solution of |1.1|) . is well defined in the phase space Hi, and, due 
to estimate (p.i|) . it is dissipative in this space. Therefore, we may speak about the associated 
global attractor. We start with the reminding of the definition of the so-calledipcally compact 
attractor which is natural for dissipative systems in unbounded domains, sec |T5f and references 
therein for the details. 

Definition 6.2. A set A C Hi, is a (locally compact) global attractor for the solution semigroup 
S{t) iff 

1) The set A is bounded in H^, and is compact in Hi^c '■= {uq G [-L^^^(R^)]^, divuo = 0}; 

2) It is strictly invariant: S{t)A = A for all t > 0; 

3) It attracts the images of bounded (in Hi,) sets in the topology of Hioc- Namely, for any 
bounded subset B C Hf, and any neighborhood 0{A) of the attractor A in the iJjoc-topology, 
there exists T = T{B, O) such that 

S{t)B C 0{A) 
for all t > T. 



The following corollary gives the existence of such an attractor. 

Corollary 6.3. Let the assumptions of Theorem life. Jl hold. Then^th 
group possesses a global attractor A (in the sense of Definition |E 
bounded solutions of (lli.l|) defined for all t G M.: 



associated solution semi- 



which is generated by all 



(6.13) 

where /C C L°°( 



A = 1C\ 



t=o\ 



jn:i 
, Hh) is a set of all solutions of ()ll.ill defined for all t € 



^ and bounded. 
WS9 



Proof. According to the abstract attractor's existence theorem, see e.g. I07"we need to verify 
two properties: 1) the existence of an absorbing set B which is bounded in Hf, and is compact 
in Hioc and 2) that the operators S{t) are continuous in the Hioc-topology on B for any fixed t. 
Indeed, due to (jb.l|l . the ball Br in the space H), will be an absorbing set for the seraiffroup 
S{t) if R is large enough although it is not compact in Hio^. However, combining (jl6.1|l with 
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, boundsattr 



s6 



6.1 



Th6 . main 



6. poly (7.2) 



6 .max 



the smoothing property (jb.3|) . we see that every solution started from uq E Br will be bounded 
in the space H^ := Hb n W^^'^iM.'^) if t > T. Thus, the i^-bah B}^ of H^ will be the desired 
absorbing set for the semigroup S{t) if R is large enough (obviously, this set is compact in the 
local topology of Hioc)- Thus, the first property holds. 

The second property is an immediate corollary of estimate (|l4.iiJ|) and the elementary fact 
that the topologies induced on B\, by the embeddings to Lf^^{M?) and Lg (M^) coincide. 

Thus, all assumptions of the abstract attractor's existence theorem are verified and, therefore, 
A exists. The formula (|b.l3|) also follows from this theorem and the corollary is proved. D 



Remark 6.4. It is not difficult to see that the factual smoothness of the attractor is restricted 
by the smoothness of the external forces g only. In articular, if g 
will be also C°°-smooth. It also worth to mention that Theorem 
that 



g|R2), then th 
land estimate 



ttractor 
snow 



(6.14) \\A\\li < C\\g\\q [\\g\\q + \\ curl^lU^^ 

where the constant C is independent of the choice of g (but, of course, depends on a > 0). 

Remark 6.5. We note that, in general, the attractor A is not compact in the initial topology of 
Hh (but only in the local topology of -f^^oc l.„JJ2^ever, if the external forces g decay as |a;| — )• oo 
{g G L'^{Q)), then arguing analogously to fTI, one can prove that A is not only compact in Hh, 
but also has the finite fractal dimension in this space. We return to the more detailed study of 
this case in the forthcoming paper. 

7. Classical Navier-Stokes problem: polynomial growth of infinite-energy 

solutions 

In this concluding section, we apply the technique developed above to the study of spatially 
non-decaying solutions of the classical Navier-Stokes problem in E?: 



(7.1) 



dtU + (u, Va 

divu = 0, ' 



lt=0 



Axu + VxP + g, 

= Uq 



which corresponds to the choice of a = in ()li.i|) . In contrast to the case a > 0, we cannot expect 
that every solution is globally bounded in time since for the simplest spatially homogeneous case 
g = const 7^ 0, we have linearly growing in time solution u{t) = tg. The purpose of this section 
is to prove that all solutions of this problem starting with uq £ Hb grow at most polynomially 
in time. Namely, the following result holds. 

Theorem 7.1. Let the assumptions of Theorem, \\b.l\ hold. Then, every weak solution of the 
Navier-Stokes problem (jlY.ip satisfies the estimate 

Mt)\\Li < Q{\\uq\\li + Mli + II cni\g\\L^){t + l)^ 



ie 



i^+) 



where the monotone function Q is independent of t and u. 



Proof. As in the proof of Theorem lb. b I we mav assume without loss of generality that curl uq S 
L°°(]R^) and we may use the maximum principle for the vorticity equation. However, the absence 
of the dissipative term aco does not allow to obtain the dissipative or even globally bounded in 
time estimate and, instead of (ib.iuli. we have 



tG 



(7.3) ||a;(t)||Loo < II curl uo||l°o + t|| curl 5^11^00, t t jr_|_, 

where the right-hand side gx.ows linearly in time. 

Similar to Theorem Id. II the parameter R in weighted energy estimates will depend on t, but 
now R{t) will grow in time in order to compensate the absence of the dissipative term au. So, 

on the time interval t € [0,r]. Following the general 



we fix a big T and consider equation 
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scheme described above, we multiply it on u^pR^xo{x) where the parameter i? will depend not 
only on the initial data, but also on T and xq G M^. Then, analogously to ()|4.3|) . we get 

1 d 



HiE (7.4) -,YM')\h..+\\^^''\\UB-)^ 



VR.x 



■^-R.^O 



< CT\\g\\l, + C{R-^ + Tf^)\\u\\l,^^j,,R^ + Ci?-i||n||i,(^,«) + \{y .P{u),m,x,u) 



where the constant C is independent on T and R. Here we have estimated the term with the 
external forces as follows: 



\{a,ULpR^x,,)\ < \\g\ 



Ll 



VR,, 



Ali < ^ll^l 



^R. 



L?„ 



+ r-i||nii2 



fR,: 



'PR,XQ 



which looks as an optimal one (in the absence of the dissipative term au) if we consider the 
solution on the time interval t G [0, T] only. 

The estimates for .the pressure term are identical to the damped case, considered before, so. 



3s tor.t 



we end up with 



6 .main 



6.mainl 



analogously to 

(7.5) ±Zn,.^{u{t))+ 

+ (2/3 - KR-^/\R-^\\u\\l2 ^ + II^IIl-)) / 0R,,,ix)\\V,u\\l2(sRs dx < CTZR^,,{g)+ 

where the constants C and K are independent on R and T . Assuming that R is large enough, 
we transform it to 

(7.6) Az^_^,^(^(t))+ 

+ (2/3 -i^i?-^/2(i?-i 11^11^2 +11,^11^^ + 1)) /" Qr^x^(x)\\V xuf^^,j^R.dx <CTZn^x,{g)^ 

\ b,R J Jx&? ^ '"" ' 

+ (c^ + KR-^I^{R-^\\u\\l2^ + M\l^ + l)\ Zr^xo{u), 

where C and K are independent of R and T. In order to be able to control the last term in the 
right-hand side, we need to assume that 



6. bad 



(7.7) 

JS.jlljainl 

Then, (py.bl) reads 



kr~^''^{r~^\\u\\l2 + ||w||loo + 1) < 2/3r-\ ie[o,r]. 



6.gr (7.8) 



^z«,,„(n(t)) - ^^ZR,xMt)) < CTZR^xoig), t G [o,r] 



and the Gronwall's inequality together with (Il2.22p gives 

6^ (7.9) ZR,xo{uit)<Ci{T + lfRi\\uo\\l2 + \\g\\l2), t e [0,T], 

1 b b 

and (|b.iu|l . we end up with 



6. good I (7.10) R-^\\u{t)\\L2 +||w(t)||L-+l < C2(r+l)(||no||i2 + ||curluo||L- + ||5'llL2 + l|curl5f||Loc + l) 



for all t G [0,T]. Estimate l/.iUl) imphes (11/. 7 1) if we fix 



6.R 



(7.11) 



R 



C2K 

/3 



|mo|Il2 + ||curluo||L°° + IblLa + ||curlfi(||Loo + Vf{T+ 1)^. 

h b 
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em6 . strange 



6. strange 



Then the approxi 
show that both I^Y 



.1 



tipn and continuity arguments (see the end of the proof of Theorem |[4.1|l 
nd (liy.iOl l ^re satisfied,^under this choice of the parameter R. It only 
remains to note that (p/.iup and (py.il|) imply (p .2^ and finish the proof of the theorem. D 

Remark 7.2. Although estimate (j^^Tessentially mMJrovesthe super-exponential upper bounds 
for the growth of u{t) in time known before (see \^\ and l[TT]j, it is probably still not optimal. 
Indeed, to the best of our knowledge there are no examples where u{t) grow faster than linear. 
However, the proved theorem gives a bit more information, namely, it factually shows that the 
mean value of u{t) over the large ball of radius R = {t + 1)^ grows not faster than linear: 

(7.12) (t + l)-4||n||i2 <C{t + l), 

b,(t+l)4 

where C is independent of time, and this estimate is already optimal since this quantity grows 
exactly linearly in time for the spatially homogeneous solution of (jly.ip mentioned at the begin- 
ning of the section. Roughly speaking, this shows that a solution with the super-linear growth 
in time (if it exists) should be "essentially non- homogeneous" in space. 



sA 



LemA. stream 



A. est 



A . path 



A . bound 



8. Appendix 1: Approximations of divergence free vector fields 

The aim of this Appendix is to construct the sequence of divergence free vector fields with 
finite supports which approximates a given divergence free vector field u G L^(R^) or u G L^(M^). 
To this end, we will essentially use the stream function G G ^/Jc ^^) which generates the vector 
field u via u = V 0. However, for general divergent free u G -L^(R^), the associated stream 



function does not belong to W^' (M ) and may grow when |x 
shows that this growth is at most linear. 

p2 



oo. The next standard lemma 



Lemma 8.1. Let u G Lf^^ 



function G W^J^ ' 



he a divergent free vector field. Then, the associated stream 
can be chosen in such way that 



.1 



|0| 






O^C I hllL2(B2^ „^) ds + C j hllL2(B2^^_^,) ds, 



'"■2) JO 

where the constant C is independent on u and (ri,r2) G M^. 

Proof. Indeed, for smooth u, the stream function G can be restored via the curvilinear integral 



^.2) 



G(x,y) 



i{x,y) 



-U2{x,y) dx + ui{x,y) dy, 



where j{x,y) = 'y{x,y,XQ,yQ) is a curve connecting an arbitrary point {xo,yQ) with (x,y), for 
instance, one may take the piece-wise linear path connecting first (2:0,^0) with (x,yo) and then 
(x,yo) with {x,y). Here (xo,yo) is an arbitrary point (e.g., (xo,yo) = (0^0)). In the general 
case when u G Lf^^{M.'^), the path integral (p. 2]) mav be ill-posed for some exceptional values of 
{xo,yo) 



so we average it over (xo,yo) ^ -Bg and write 



A.pathl (8.3) Q{x,y) 



1 



m 



{xo,yo)<^Bl \J'y{x,y,xo,yo) 



-U2{x,y)dx + ui{x,y)dy] dxodyo. 



It is not difficult to check that this integral is well-defined for any u G -L^^^(R^) and indeed 



u = V"'-G for any divergent free u. Moreover 
and the lemma is proved. 



T|) also follows by the straightforward estimates 

D 



Corollary 8.2. Let u G L^(M^) he divergent free. Then there exists a stream function G G 
Wl^l{M?) such that 

(8.4) 



(|xo| + 1) ^||G||i2(Bi^) <C\\u\\l2, 
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A. dot 



CorA . app 



app 



A.unlim 



sA2 



efA2. dyadic 



A2.PL 



A2.0 



A2.HPL 



where C is independent of xq. If, in addition, u G L'^{M.'^), then 



i.5) 



M 



lim (|xo| + l) ^||e||i2(m ) 



0. 



Indeed, both statements are immediate corollaries of estimate (jb.i 
The next corollary can be considered as a main result of the Appendix. 

Corollary 8.3. Let u G L^(M^) he divergence free. Then, there exists a sequence u 
N gN, of divergence free vector fields such that 



^elK 



'), 



(8.6) 






< C\\u\ 



N 



Li^ 



suppti C i?, 



,2Af 
' 



u^{x) = u{x), xeB^, 



where the constant C is independent of N and u. If, in addition, u £ L^(M^) then 
(8.7) 



lim 

|xo|^>oo 



1^ WlHbu 



uniformly with respect to N . 
Proof. We may fix 



u^^(x) := V^(e(x)c^7v,o(a;)) = 'u(x)(/^ij„o(x) + e{x)V^ipR^o{x) 



!to e(x) iyh 

(15^ and i^m. 



^_^-. stream 

e stream function constructed in Lemma 118.11 and (fn 



,xoy 



are defined via 

J. Indeed, since (pN,o{^) = 1 for \x\ < N, we have u^{x) = u{x) for such x 

and ^pNfi{x) = for \x\ > 2N implies that u^{x) = for such x. Moreover, due to (|l3.16p . 
^^Ro(x)\ < CN^^ and this gradient is non-zero only if A^ < Ixl < 2N. Thus, the uniform 

' lanri lA .hm mH ' ' ~n\ .im 1 im ^„ „, 

in the L^(]S 1 follow from 



bounds 



and the uniform limit 



fa2£ 



1.1 1 II 1111 • I-.. 

anfe"r u G Lli 



immediate corollary of ip.b\] . Corollarv ll8.ai is proved. 

9. Appendix 2: The interpolation inequality 



IS an 

D 



The aim of this Appendix is to verify the sharp interpolation inequality stated in Lemma Ib.bl 
Although this inequality looks more or less standard, it is not easy to find the precise reference 
in the literature, so for the convenience of the reader, we sketch its proof here. To this end, 
we need to remind briefly the definitions arul some facts from the theory of Besov spaces and 



Paley-Littlewood decomposition, see l[13fl2U] for more detailed exposition. 

Definition 9.1. Let G C^(M ) be a non-negative cut-off function such that (/)(^) = 1 for 
1^1 < 1/2 and ^(C) = for |^| > 1 and let ^(6 = </'(C/2) - 0(0- Then, for every j G Z, we 



4 . key 



define operators Sj, A,- : S' 



S' 



as follows: 



S,f:=<PiC/2^)f, A,f:=iP{C/2^)f, / G 5'(M'^), 

where / is the Fourier transform of the tempered distribution / G 5'(M ). Then, for any N £ Z 
andany /g5'(M'='), 

(9.1) f = SNf+Y.A,f. 

The distribution A^ is called jth dyadic block of the distribution / in the Paley-Littlewood 
decomposition l^.i\] . If. in addition. 



(9.2) 



N 



lim SNf = in S'{R'^) 



then the homogeneous Paley-Littlewood decomposition holds: 
(9.3) f = ^^jf- 
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A2.Ber 



DefA2.Besov 



A2.Bes 



A2.HBes 



The space of distributions satisfying 
and is denoted by S^^iR'^). 



is called the space of distributions vanishing at infinity 



PL 



Obviously, L'^{M.'^) C 5g(M'^) and by this reason the homogeneous decomposition (jlu.3|l holds 
for all functions used in this Appendix. We also remind that the operators Aj commute with 
differentiation and satisfy the Bernstein type inequalities: 



(9.4) 



1) \\Ajf\\LP<C2-mVAjf\\LP, 

ltS2, 
where I < p < q < oo and C is independent of j, see e.g., 



Definition 9.2. For any 1 < p,q < oo and any a G IR+, we define the Besov space Bp {M. ) as 
the subspace of iS'(M'^) generated by the following norm: 



1/9 



(9.5) 



B,„ := llSifhr' + \^i2'''\\A,fh,) 
vi=i 



< oo. 



The homogeneous Besov space Bp^g(R'^) is defined via 
(9.6) 



i/g 



B" 



J2 (2'^^||A,/||,^^'^ 



LP) 



< OO. 



PropA2.ZC 



A2.Riesz 



^ IS only a semi- norm and the space B^^ is defined by modulo of polynomials 

'o( 



In contrast to .i,^.^,, .^ ^.^., ^ ^^^... x.v^^x.^ ^^.^ ^.^^ ^^^^^ ^„„ 
|L02 ^ ^ -^ ^ p,g 

(see iTJj for the details). However, it is not essential for us since it is a norm on Sq(M ) and we 



will use these spaces in the situation when / G L (R ) C Sq(M ) only. 



Finally, the space BMO 
following semi-norm: 



of functions with bounded mean oscillation is defined by the 



(9.7) 



BMO 



•= s^P T^ / \f{x)-mBf\dx < CO, niBf ■■= 7-^ f{x)dx, 
BeB \B\ Jb \B\ Jb 



where ;B is a collection of all balls B^ , < i? < oo, xq G M and |i?| is the d-dimensional volume 



of the ball B. This space is again defined by modulo of constants and it is again not essential 
for us since we will work with functions from L^(M ) Pi BMO(^ ) only. 

The next standard proposition is crucial for what follows. 

Proposition 9.3. Let Rj, j = 1,- ■ ■ ,d, be the Riesz operators defined by 



(9.8) 



R,f 






/• 



to 



Then the compositions of Riesz operators Rij := Ri o Rj are bounded operators from L^ 
LP(M'^) and from L°°{W^) to BMO{W^). Moreover, 

BMO(M'') C i?So,oo(K'') 

and, consequently, Rij are bounded operators from L°°(]R ) to B^ ^{M. ) as well. 

|r.02 I LailL4 . key 

For the proof of this result see, e.g., \T3. We are now ready to prove Lemma lb. 6 1 

|r.em 4.kev M^Jjayest J 4^Js£yestl 

Proof of Lemma \\d.t)[ We first note that (|b.i3|) and (jb.i4|) are scaling invaria|nt, so we only need 



Jr ■ I .^7' 

to check them for, say, R = 1/2. We start with the most delicate estimate ()b.lci|l 
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Extending the vector field by zero outside of Bq, denoting hi = divn and /i2 = curlu, after 
the Fourier transform, we get 



A2.Helin 



(9.9) 



6 



6 



^1 = ~M T9 79^1 + 79 79"2 , U2 = -t\ -n -n/ll 



e?+cr ef + er^' 



6 



^'^) 



B.2 



,2 / ^ei+e^ a 

and we see that V^^t can be expressed through hi and /i2 via the compositions Rij of Riesz 
operators. Thus, due to Proposition fe.ai V^n £ [i?J^ oo(M^)]^ and 

(9.10) IIVxull^o^ ^-^ (||^iIIl°°(b1) + II^2||loo(bi)J . 

Thus, using the Bernstein inequahties 

IIAjuIIlp < C\\Aj\7xu\\lp2~\ 
we see that u G BL, „(M^) and 



B.3 



(9.11) 



M lii 



Bi (IR2) < C'(||/ii||loo + ||/i2||L°°) • 



At the next step, we interpolate between B22 and B^ ^. Namely, by the Holder inequality 



and, therefore. 



All ^ II A l|5/6|| A i|l/6 



5/12 



B.int 



(8-12) ll«ll,t ,.„, = E(2"»l|A,.||,„.)-'= < 



12 12 



,jez 



5/12 



^^ Y^iiA ||2 /ojiiA II A^/S \ II ||5/6 II ||l/6 



rjO ll"'llr;l 



int. 



B.besov-hom 



Note that, due to the Plancherel equality, L'^{M.'^) C B^^^O^"^)^ so, combining (||jj.'l2f) and (jl9^ . 

we have 

(9.13) iu,ii , ^ r'lu, 11^/6 ni^ju^ -L ii^„iu_^V6 






12 12 



< C||n||^'2 (||/ii||l°° + ||/i2||L°°) 



B.ZC 



.13|l for the usual (non-homogeneous) Besov 
spaces. To this end, we utilize the fact tlmt the supports of hi, /i2 and u belong to the unit ball 
Bq. By this reason, Proposition ib.^l gives 

(9.14) \\u\\l2{bI) < C'||V^'u||i2(Bi) < C'l (||/ii||l2 + ||/i2||l2) < C2 (||/ii||l- + ||/i2||L-) 

Therefore, due to the Bernstein inequalities, 

||S'l'u||j;^12/5 < C||S'lli||i2 < Cl||n||i2 < C2 (||/il||Loo + ||/i2||L°°) 

/ iP. . 9i . [R .hfipov-hom 

and, together with ()b.l4p and (lb.i;^ll . we have 



B. besov 



(9.15) 



\u\\ 1 

B 



1 5/6 



12 12 



< C||n||^2 (||/ii||l°° + ||/i2||L°°) 



1/6 



Using now that Zg^ C L^ for q2 > Qi, we have Sfa 12 C i?i2 and 



,3 



B.besovl 



(9.16) ||li|| 1 ^ >-^ll"ll r2 

^ ' " "b?2 3(K2)" "^ 

Finally, using the embedding theorem 



<C||^.f/^||/il||L- + ||/^2||L=o)l/^ 



B.Triebel 



(9.17) 



^'^ p q 
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ITriTS J lJiayest 

(see |'2tT|7'page 206) with s = 1/6, n = 5/12 and q = 3, we end up with (jb.l3|) . 

Thus, it only remains to verify (jbJjlL m the particular case xq = and 2R = 1. Indeed, due 
to Proposition lb.;!?! and expressions (|y.9[) . we have 



ll^llwi>p(Bi) ^ C\\Vxu\\^p(^Bi^ < Cidl div u\\lp + II curlnllip) 
and estimate ()b.i4lj is now an immediate corollary of the interpolation inequality 

II^iIloo < c||'u||^2||u||J;^;i,p, 

|Trj^78 I LaiIL4 . key 

see pCTfTThus Lemma lb.6l is proved. D 
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